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Abstract 

The change of the probability of proton emission in nuclear photoeffect due to an intense coherent 
(laser) field is discussed near the threshold, where the hindering effect of the Coulomb field of the 
remainder nucleus is essential. The ratio of laser-assisted and laser free differential cross section is 
deduced and found to be independent of the polarization state of the 7 field and the two types of 
initial nuclear state considered. The numerical values of this ratio are given at some characteristic 
parameters of the intense field. 
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I. INTRODUCTION 



The development of constructing coherent electromagnetic sources of higher and higher 
intensity with increasing photon energy up to the hard x-ray range motivates the theoretical 
study of the change of the processes of strong.y bound systems, such as the nuclear processes, 
by these intense fields In this paper the change of the nuclear photoeffect due to the 
presence of an intense coherent electromagnetic field is studied. This process is analogous to 

;he laser-assisted x-ray photo effect (x-ray absorption), a process which was already discussed 

2| in the late 80 's taking into account gauge invariance 

a, a 

. The laser-assisted nuclear 

photoeffect (LANP) and the laser-assisted x-ray photo effect (x-ray absorption) are laser- 
assisted bound-free transitions. The difference between them lies in the charged particle 
(proton or electron, respectively) which takes part in these processes. Although the LANP 
was recently investigated far from the threshold and neglecting the effect of the Coulomb field 
of the remainder nucleus Q but in the case of the laser assisted x-ray absorption processes 
it was found that the most interesting changes due to the presence of the laser field appear 
near the threshold [6fl, Ml. 

n 

Thus, applying the results of [2[ the LANP is reexamined in a gauge invariant manner 
and near the threshold, where the hindering effect of the Coulomb field of the remainder 
nucleus is very large so that it must be taken into account. The effect of the Coulomb field 
of the remainder nucleus on the transition rate is approximately taken into account. The 
laser-modified differential cross section is compared to the laser-free differential cross section, 
and it is shown that their ratio does not depend on nuclear parameters in the two types of 
initial nuclear states investigated and on the state of polarization of the 7 radiation but it 
has only laser parameter dependence. 

The process investigated can be symbolically written as 



10. 



7 + nco + illY -> A Z X + ]p, (1) 

where z+]X denotes the target nucleus of mass number A + 1 and of charge number Z + 1. 
The target nucleus absorbs a 7 photon symbolized by u 1 and n number of laser photons 
take part in the process which is symbolized by uuq. n < and n > correspond to |n| 
number of laser photon emission and absorption, respectively. As a result a free proton \p 
is emitted and the remainder nucleus is %X. 
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The calculation is made in the radiation (pA) gauge and in the long wavelength ap- 
proximation (LWA) of the electromagnetic fields, the recoil of the remainder nucleus and 
the initial momentum carried by the laser and 7 fields are neglected. In the case of a 
circularly polarized monochromatic wave for the vector potential of laser field Ai(t) = 
A [cos((jj t) e 1 — sin(o;ot) e 2] is used. oo is the angular frequency of the laser. The ampli- 
tude of the corresponding electric field E = ujqAq/c. The frame of reference is spanned by 
the unit vectors e x = e 1, = ~£ \ and ~~ef z = ~e\ x ~£ 2 - The vector potential describing 
the gamma radiation is exp (—iu-ft), with hu-y the energy and 1^ the 

unit vector of state of polarization of the 7 photon, and V the volume of normalization. 



II. GAUGE INVARIANT S-MATRIX ELEMENT OF LASER-MODIFIED PRO- 
TON EMISSION IN NUCLEAR PHOTOEFFECT 

It is shown in (^J that the electromagnetic transition amplitudes of a particle (proton) of 
rest mass m and of charge e in the presence of a laser field are determined by the matrix 
elements of the operator — with the eigenstates of the instantaneous energy operator 

e9=^(f-~Atf + V(r) (2) 

in both (rE and pA) gauges, (e is the elementary charge and the superscript g refers to the 
gauge.) Accordingly the gauge independent S-matrix element can be written as 

i 



Sfi ~~h 



JdtJ d 3 rr f (-e^ • t (t)) ^, (3) 



where ipi and ipf are the initial and final states of the proton in the same gauge and h is the 
reduced Planck constant. 

Our calculation is carried out in the radiation (pA) gauge because of the choise of the 
final state of the proton (see below). The initial state of the proton has the form 

^ = e(^)^)e-^*, (4) 

where 0o(~^) is a stationary nuclear state of separation energy of the proton. The 
exp (i e ^, ) factor, where ~~X = A*l (t) + appears because of gauge transformation 
since 0o is the eigenfunction of the instantaneous energy operator 

e E = ^f 2 + V N {r) + V iC {r) (5) 



in the rE gauge. Vjv(r) is the nuclear potential and Vic(r) is the Coulomb potential felt by 

the proton initially and the superscript E refers to the rE gauge. The modification of the 

initial state due to the laser field is neglected since the direct effect of the intense laser field 

fl 

on the nucleus has been found negligible [81] at the laser parameters discussed. It is also 
supposed that the initial nucleus does not have excited state which is resonant or nearly 
resonant with the applied 7 radiation. 

If similarly to [5j the modification of the final state due to the strong interaction is 
neglected then in the final state and in the pA gauge the instantaneous energy operator e R 
reads as 

e R =^(t-- c Al(t)) 2 + V c (r), (6) 

where the superscript R refers to the radiation (pA) gauge and Vc(r) is the Coulomb po- 
tential of the remainder nucleus. 

An approximated solution of (jSJ), i.e. an approximated time dependent state of a particle 
in the laserplus Coulomb fields is the Coulomb- Volkov solution of a proton of wave number 
vector (J (g|, |lo|: 

^ (?,t) = l^V^x^^exp (-iEt/H) f(t). (7) 

Here V A_1//2 e 1 ^ '^x(^ ^) is the Coulomb function, i.e. the wave function of a free proton in 
a repulsive Coulomb field of charge number Z, V denotes the volume of normalization, r 
is the relative coordinate of the two particles. 

X($,?) = e-^ 2 T(l + irj)iF x (-i<n, 1; i[Qr - ~$ ■ !>]), (8) 

where 

n(Q) = z a ,™ (9) 

is the Sommerfeld parameter, with a / the fine structure constant and it is supposed that m 
much less than the rest mass of the remainder nucleus. 1F1 is the confluent hypergeometric 



function and F is the Gamma function 11]. 
The function 

fit) = exp[ia sin(w t + 770)] (10) 

where 

a = a$ sin (1?) with a$ = — (11) 



mu. 







Here the polar angles of the wave number vector of the outgoing proton are •& and 770, i.e. 
they are the polar angles of the direction in which the proton is ejected. 

In the low-energy range (QR <C 1, where R is the radius of a nucleon) and for \~r^\ < R 
the long wavelength approximation yields 



X0^)^ = Xc(Q) = J =- , (12) 

7^=0 y exp [2717] (Q)\ - 1 

which is the square root of the so-called Coulomb factor. (The Coulomb factor (y 2 c {Q)) 
describes well e.g. the Coulomb correction to the spectrum shape of beta decay |12|.) 

For the final state of a proton of wave number vector the LWA of the nonrelativistic 
Coulomb- Volkov solution ip-^ is used. 

^,lwa 0*, t) = Xc(Q)V~ 1/2 e^^ exp (-iEt/h) fit) (13) 

with E = h 2 Q 2 1 (2m) + U p , that is the energy of the outgoing proton in the intense field, 
where U p = e 2 Eq / '(2mujQ) is the ponderomotive energy. 

Substituting (jl]), (fT3"j) into ([3]) and using one can obtain the following form 

of the gauge independent S-matrix element as 

S fl = J eX p[* (E + E b ) t/h}f*(t)^G [f (t)} dt, (14) 

where 

G(t) = J U^^^dh (15) 
is the Fourier transform of the initial stationary nuclear state 0o( ' ) of the proton and 

t(t) = cl--^2. (16) 

(Equation (fT4|) can be obtained directly with the aid of Eq.(27) of j^J.) 
Using the d t G = (d q G) £ (d A .q) (d t Aj) identity 

d ^ / d „\ e -3 e 



dt 



with = -\d t ~t, i.e. if = ~t L (t) + 

The Q ■ E L term of the last factor of (1TT1) can be neglected if the pure intense field 
induced proton stripping process is negligible since this term describes the process without 
the gamma photon. Furthermore, the ratio of the amplitudes of A* y ■ E L and A l • ^ 7 equals 



wl/u; 7 <C 1. Therefore the (jf • E — j-~k-E* = ■ i^ 7 — j-^^ - E*^ approximation is justified 
to use, where i^ 7 = i^2itFu J j 1 /Vlt exp (—iu^t). The relative strength of the 7 and 

7 terms is characterized by the parameter 5 = eAo/ (hcQ). In the laser free case 
Q = \ / 2mA/h, where A = hu y — E b is the difference of the photon energy and the proton 
separation energy. Numerical estimation shows that 5 ~ 0.05 near A = 50 keV used here 
and in the case of laser photon energy and intensity values discussed. Therefore the 
term is the leading one in the last factor of ( ITT)) . As to the radiation field dependence of 
d q G 1 the effect of y£ 7 is negligible in ~~<t(t) and thus ~q* (t) = Q — wzAl- It was shown 
above that the amplitude of oscillation of q (t) due to the intense field can be neglected. 
Therefore q = Q can be used in d q G, and ( 1181) results 



Using the Jacobi- Anger formula in the Fourier series expansion of /*(£) [13| the S-matrix 
element can be written as 

S fl = f 2 ^f )] \ c(Q) x (19) 

n=no 

where £ =Q-1?/Q, J n (a) is a Bessel function of the first kind, and 

m\ hQ 2 , u p + E b 

u n \.Q) = ^ — H ^ w 7 - nw . (20) 

The terms, that are small if the <C 1 condition is fulfilled and /3 2 see below), 

were neglected in the calculation. 

III. GAUGE INVARIANT DIFFERENTIAL CROSS SECTION OF LASER- 
MODIFIED PROTON EMISSION IN NUCLEAR PHOTOEFFECT 

In this paper two cases of initial nuclear energy — E b of the initial state having different 
type of space dependent part 0o (~^) are considered in order to show the general nature of the 
effect of the laser on the process. The one case is the 8 B one-proton halo isotope of separation 
energy E b = 0.137 MeV Q and of initial state O (^) = (2tt) _1 p\ ,2 e~ Pir / {for), with 



6 



0i = fi^2mEb^/h and m the rest mass of the proton (/i = 1.84, 0i = 1.495 x 10 12 cm 1 ) 



In the other case the initial state is the vs\/2 shell model state [l5[| of the form 0o 



N v e~?l 2 F (— i/,!,/? 2 ), with iV„ = /3 2 3/2 [r(z/+ |)/(27rz/!)] 1/2 /r(|) where ^ = 0,1,2,... is 
the quantum number of the nuclear shell model and T(x) denotes the gamma function. 
F (—u, |,p 2 ) is the confluent hypergeometric function, p = /3 2 r, /3 2 = ^Jmuish/K m is the 
nucleon rest mass and uj£h is the shell model angular frequency (hush = 40v4~ 1//3 MeV, A 



is the nucleon number 



To), and (3 2 = 9.82 x lO 12 ^- 1 / 6 



cm 



The differential cross section of LANP has the form 



da 
dVLn 



oo , 



dQ„ 



(21) 



n=no 



where dfl q is the differential solid angle around the direction of the outgoing proton. n < 
is the smallest integer fulfilling the A + nftwQ — U p > condition and A is the same in both 
cases of the initial state. (The cases n < and n > correspond to |n| number laser photon 
emission and absorption, respectively.) 
The partial differential cross section 

da r . 



dQ n 



Vo,n (Qn) |£| JjatfnSin^) 



with 



CQ.n (Qr 



(Q n ) [d Q G 



(22) 



(23) 



Here Q n = ^J2m [A + nftwo — U p }/h, £ = Q n -~^ ' /Qn, ^ P is the reduced Compton wavelength 
of the proton and fc 7 = co 7 /c. 



Q, 



X, 



2A 



mr 



with e n 



nfroj — U„ 



1 + 



A 



(24) 



J n (c(^ n sin {}) is a Bessel function of the first kind, with a$ n = eE Q n / (touqJ 
In the variable e n the Coulomb factor Xc rea ds as 

K Cb 



Xc(£n) 



exp 



K, 



1 



(25) 



where Kcb = 2ti Z a j \f mc 2 / (2A) with Z the charge number of the remainder nucleus. The 
Coulomb factor causes a strong hindering of the effect in both the laser assisted and laser 
free cases. 
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FIG. 1: Laser photon number dependence of R n (see (|30p ). that is the rate of change in one channel 
due to the presence of intense field. The cases n < and n > correspond to |n| number laser 
photon emission and absorption, respectively. The charge number of the remainder nucleus Z = 4, 
A = hu^ — Ef, is the difference of the gamma photon energy and the proton separation energy, 
i? = vr/2, the laser intensity I = 10 20 Wcm~ 2 and the laser photon energy fkjQ = 1.65 eV. 




50 60 
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FIG. 2: The A dependence of R (see (|2U]) ) at i9 = vr/2 in the case of Z = 4 with J = 10 20 WcrrT 2 
and fajn = 1-65 eV . A = /io; 7 — Eb is the difference of the gamma photon energy and the proton 
separation energy. 

Near the threshold (Q n (3±, fa) the [OqG]^ oc e 2 n in the case of the two types of initial 
state discussed and 

An 

(26) 



dVt n 



S 
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FIG. 3: The i? dependence of R n (see (J3QK) with I = 10 20 Ifcm" 2 and Hlu = 1.65 at (a) 
n = 18000 (b) n = 14000 (c) n = 10000 and (d) n = 5000. The charge number of the remainder 
nucleus Z = 4 and A = 50 keV. 




0.5 1 1.5 2 2.5 3 
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FIG. 4: The dependence of i? n (see fl30])) with n = 10000 and ^ = 1-65 eV at (a) I = 4 • 10 19 
Wcm~ 2 (b) 7 = 6- 10 19 Wcm~ 2 and (c) I = 10 20 Wcm~ 2 . The charge number of the remainder 
nucleus Z = 4 and A = 50 fceV. 

Here o"o is constant, it depends on the form of the initial state and 

e 2 n J 2 (ay n sin??) 
exp (K Cb /e n ) - 1 

In the laser free case (a$ n = 0) using e (kp = 0) = 1, J, 2 (0) = 1 and J 2 (0) = at n ^ 0, 
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the differential cross section near above the threshold 



da th 
dVL q 



aoSth Kl 2 (28) 



with S th = [exp (K C b) ~ 1] i - 

IV. NUMERICAL RESULTS 

The ratio R of the laser-assisted and the laser-free differential cross sections, 

oo 

R=Y, R n (29) 

n=no 

where R n = S n / S t h [ (l26|) divided by fl28|) ]. equals the ratio of the rates of the corresponding 
processes in an elementary solid angel in a given direction of the outgoing proton. The rate 
of change in one channel 

exp (K Cb ) - 1 2 j2 . . . . 

Rn = 7T7 / v T £ n J n i a $n Sin it) (30) 

exp (Kcb/£n) - 1 

with 



a$ n = e n hceE (ftw ) \/2A/ (mc 2 ) (31) 

in the variable e n . So _R n and R describe the change caused by the intense coherent field 
independently of the state of 7 polarization and the initial states applied. 

In our numerical calculation the laser photon energy TujJq = 1.65 eV. First the case 
when the outgoing proton moves in the plane of polarization of the laser beam (i? = 7r/2 
(sin?? = 1)) is investigated in the case of Z = 4. Fig. 1 shows the laser photon number 
dependence of R n at A = 50 keV with laser intensity / = 10 20 Wcm~ 2 . The intensity 
dependence of R has been investigated with A = 50 keV. R linearly increases from R = 1.00 
at / = 10 18 Wcm~ 2 up to R = 1.12 at / = 10 20 WcmT 2 . Fig. 2 depicts the A dependence of 
R with I = 10 20 Wcm~ 2 . Fig. 3 shows the d dependence of R n in the case / = 10 20 WcmT 2 
at (a) n = 18000 (b) n = 14000 (c) n = 10000 and (d) n = 5000. Finally Fig. 4 depicts the 
$ dependence of R n at n = 10000 with different laser intensities: (a) / = 4 x 10 19 WcmT 2 
(b) I = 6 x 10 19 WcmT 2 (c) I = 10 20 Wcm~ 2 . The numerical calculation in all the cases 
discussed above has been repeated at Z = 9 and negligible change has been found. 
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V. DISCUSSION AND SUMMARY 



In order to compare our results and the results of 5| we investigate our formulae in the 
hujj 3> E b limit. In this case e n — > 1, A — > hui 7 , 



a$ n — > oc# — hceE (Huj ) y2huj^/ (mc 2 ) (32) 

and 

Rn ^ J n 

(a#sin#). (33) 

Thus in this limit the n dependence disappears from the argument of the Bessel function 
and 

oo 

R^Yl ^nKsin^) - !■ ( 34 ) 

n=n 

The cross sections obtained in |5| are symmetric in n around n = (see the Figures 
of j^j) and the total cross section is asserted to be uneffected by the laser radiation. This 
corresponds to R — 1. In contrast, our result is significantly asymmetric in n ( see Fig. 
1.) which is the consequence of the e n dependence of R n (see (|30|) ). The change (increase) 
of the kinetic energy of the proton is manifested in the increase of e n from e_nooo = 0.761 
up to e isooo — 1-240. The increase of e n with increasing n causes an asymmetry in the n 
dependence of Fig. 1. The sum of the changes (increments) in the different channels results 
a moderate increment of R (R < 1.28) as can be seen in Fig. 2. 

Summarizing, one can say that near the threshold R, which measures the change of the 
rate of the LANP, has minor A and intensity dependence, and negligible Z dependence. 
Furthermore, R n , which is the rate of change in one channel (at a definite laser photon 
number), has significant laser photon number and $ dependences. R and R n are the same in 
the cases of the two different initial states considered. Since the numerical results obtained 
seem to be independent of the initial nuclear states chosen it can be expected that R and 
R n have minor dependence from the form of the initial state in general. 

Regarding real experimental situation at such a high intensity supposed it is hard to 
distinguish photo-protons from background protons that arise as a result of interaction with 
the hot electron plasma created by the intense laser field. We also have to mention that in an 
experiment the gamma ray pulse from an accelerator must be synchronized with an intense 
(e.g. attosecond) pulse of a laser system. The wispy target determined by the focal spot of 
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the focused intense laser beam, the low repetition rate of the laser system, and the angular 
resolution of the proton detector together make huge challenge to carry out a successful 
near threshold, laser-modified proton emission experiment that could produce significant 
counting statistics. 
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